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Abstract 

We present a study of entanglement in the case of the ID extended anisotropic Heisenberg model. We investigate two quantum 
phase transitions (QPTs) within the previously found ergodicity phase diagram [E. Plekhanov, A. Avella, and F. Mancini Phys. 
— — ,Rev. B 74, 115120 (2006)]. Our calculations are done by means of the numerically exact Lanczos method at T — 0, followed by a 
^ finite-size scaling. As a measure of entanglement we use the concurrence for two spins out of the system. We conclude from our 
d studies that these QPTs are accompanied by a qualitative entanglement change. 



^\ Key words: Entanglement, ID Spin systems, Heisenberg model, Lanczos 
rH ~PACS: 75.10.Dg, 75.10.Pq, 75.30.-m, 75.30.Kz, 75.50.Ee, 77.84.Bw 
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1. Introduction 

Quantum phase transitions (QPTs) are among the most 
fascinating topics in quantum mechanics. Recently, in con- 
nection to the quantum information theory (QIT), the cor- 
relation between QPTs and entanglement has been exten- 
sively studied [1]. Many systems have been found where 
QPT is accompanied by a qualitative change of the en- 
tanglement, suggesting an implicit connection between the 
two. The spin systems are a natural implementation (at 
least theoretically) for QIT devices due to the isomorphism 
between the Hilbert spaces of a single spin and the one of a 
qubit - the central object of quantum information. From the 
QIT point of view, the spin chains represent a non-trivial 
example of the qubit network, and it is therefore legitimate 
to look for the entangled states of such systems. There ex- 
ist a few entanglement measures, which differ mainly in the 
way the system is split into blocks, whose entanglement is 
measured. The one-tangle, or von Neumann entropy, is only 
a function of the local magnetisation and, hence, is scarcely 
informative. It is therefore more preferable to use the con- 
currence [2], or pairwise entanglement, which, on the con- 
trary, depends on the spin-spin correlation functions. The 
concurrence for a couple of spins at sites i and j is defined 
as: 



Cij = max(0, Ai — A2 — A3 — A4), 



(1) 
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where {A,} are the eigenvalues, in decreasing order, of the 
Hcrmitian matrix R, which can be expressed in terms of 
the two-point correlation functions [2]. Its extreme values, 
zero or one, indicate that the system is either a product 
state or maximally entangled, respectively. In order to take 
into account the contribution to the entanglement, coming 
from the correlations at all distances we adopt Ti [3] as a 
measure of the strength of the entanglement: 
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(2) 



where Ci^+d is the concurrence between the site i and its 
d-th neighbor. In a translationally invariant state T2 does 
not depend on i. 

2. Model and Method 

We measure ti in the one-dimensional anisotropic ex- 
tended Heisenberg model with next-nearest-neighbor inter- 
action on a chain with L sites, subject to periodic boundary 
conditions: 



H — —J z 3fSf_ 
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We diagonalize the Hamiltonian (3) by means of the Lanc- 
zos diagonalization technique. In doing that we take into 
account translational symmetry and classify the eigenstates 
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(a) Phase diagram of (3) in the J' — J± plane. 
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(b) Transition A: J' is set to 0.2J Z 
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(c) Transition B: Jx is set to 0.5J Z 



Fig. 1. (a) Schematic phase diagram of the model (3). See Ref. [4] for details. Notice the locations and directions of the lines crossing the 
two phase transitions. Entanglement measure T2 across transition A (a) and transition B (b) for different system sizes L = 12, 16, 20, 24. 



by the eigenvalues of total S z which is a good quantum 
number. We use J z > 0, which corresponds to ferromag- 
netic coupling. 

3. Results 

The phase diagram in the J' — J± plane, as found in 
Ref. [4] and shown on Fig. 1(a), manifests the presence of 
a totally polarized phase (NE-II phase), which neighbors 
two paramagnetic phases at the lines J± = J z and J' « 
0.33 J z . In this article we study these two phase transitions, 
which we will call transition A and transition B, respec- 
tively. In the NE-II the ground state is doubly degenerate 
with all spins either "up" or "down" and therefore r 2 = 
because this is a product state. The behavior of r 2 across 
the transition lines for two representative values of J' and 
Jj_, is shown in Fig. 1(b)- (c), respectively. We summarize 
the properties of the two transitions as follows: 

Transition A: In this transition frustration, induced by 
the antifcrromagnetic interaction, parametrized by J± , de- 
stroys the ferromagnetic order. If J± > J z the system be- 
comes antifcrromagnetic in XY plane, as emerges from our 
analysis [5]. It is clear that such a ground state is entangled. 
C^j+d monotonically decreases as a function of d, and the 
main contribution to t 2 comes from the nearest-neighbor 
term (d = 1). For the system sizes considered, L = 12 24, 
the decay of the antiferromagnetic correlations is very weak 
and they scale well with increasing L. This fact explains 
the good finite-size scaling of our results for r 2 . 

Transition B: In the definition of the concurrence (1) it 
might happen that for some d there is no dominant eigen- 
value of R and Cj^+d = 0. This would mean that for this 
particular d the two spins are untangled from the rest of 
the system. This is the case of the phase, situated behind 
the transition B. Precisely, the only non-zero contribution 
to t 2 comes from d = 2. This implies the enhancement 
of correlations at d = 2, as indeed found by examining 
the correlation functions. This feature is short-ranged even 
for the small sizes considered here and persists both in 
Z— direction and in XY plane. Within the range of J' con- 
sidered [0.3J Z ,0AJ Z ] there are regions, characterized by a 



finite magnetisation per site less than 1/2. For J' > 0.37J 2 
the magnetisation stabilizes at zero. It is this finite mag- 
netisation which is responsible for the steps of t 2 on the 
right panel of Fig. 1(c). Unfortunately our computation fa- 
cilities did not allow us to conclude whether these steps are 
a finite-size effect, and whether the vanishing of the ground 
state magnetisation at the QPT is continuous or not. The 
shorter range of the correlations respect to the Transition 
A case and the presence of the magnetisation steps imply 
a worse finite-size scaling in this case. 

In conclusion, we have studied the entanglement change 
across two QPTs in the vicinity of the ferromagnetic phase 
in the extended anisotropic Hcisenbcrg model (3). The en- 
tanglement appears to be intimately connected to the un- 
derlying ground state of the system. Starting from a totally 
untangled ferromagnetic state, by changing the Hamilto- 
nian parameters, we observed the complete reordering of 
the ground state structure. Such reordering is accompanied 
by the appearance of pairwisc entanglement, reflecting the 
structure of underlying ground state. Namely, in the tran- 
sition A, entanglement is non-zero for all the distances, be- 
ing maximal for nearest neighbors, while for transition B 
only the next-nearest-neighbor entanglement is non-zero. 
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